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Abstract: The inflationary flow equations are a frequently used method of survey- 
ing the space of inflationary models. In these applications the inflnite hierarchy of 
differential equations is truncated in a way which has been shown to be equivalent to 
restricting the set of models considered to those characterized by polynomial infla- 
ton potentials. This paper explores a different method of solving the flow equations, 
which does not truncate the hierarchy and in consequence covers a much wider class 
of models while retaining the practical usability of the standard approach. 
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1. Introduction 

The inflationary flow equations, introduced by Kinney following earlier work by 
Hoffman and Turner p| are frequently used as a means of surveying the space of 
scalar field theories describing inflation. The flow equations form an infinite hier- 
archy of ordinary differential equations. They provide a convenient framework for 
parameterizing the space of inflationary solutions of Einstein's equations coupled 
to a single inflaton with canonical kinetic terms. These equations form the basis 
of Monte-Carlo reconstruction of the inflaton potential [Q], as well as other studies 
which explore the space of inflationary models [|^-||14||. 

Practical applications of the flow equations involve truncating the infinite hier- 
archy so as to obtain a closed set of equations which can be solved. The procedure 
introduced by Kinney 0] and used by all subsequent studies^ defines subspaces of 

^Some interesting alternative approaches to characterizing the space of inflationary models can 
be found in the paper of Ramirez and Liddle |Tq| . 
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solutions characterized by all Hubble flow parameters vanishing apart from a finite 
number. It was shown by Liddle [l^ that this procedure corresponds to restricting 
the set of all infiaton potentials to polynomials of order related to the number of flow 
parameters allowed to assume non-zero values. 

The flow equations themselves do not make any assumptions about the potential 
energy density which defines a specific model. Each solution of the flow equations 
however corresponds to some definite potential, which can easily be obtained. Thus 
subspaces of the space of all solutions correspond to definite classes of potentials; 
dynamical information enters the flow equations algorithm by the means chosen for 
truncating the hierarchy. That is the point when the class of potentials to be scanned 
is determined. 

While straightforward (and adequate for many purposes), the truncation scheme 
introduced by Kinney and universally employed in subsequent studies of the flow 
equations excludes some interesting models of inflation, for example those involving 
exponentials of the infiaton field. Such cases arise in some supergravity or string 



motivated models [p7f| , pTq] , and so it would be nice to be able to broaden the scheme so 
that they could be included. One may argue that for studying physical effects which 
are sensitive only to a limited range of infiaton values a polynomial approximation 
for the potential may be all that is needed, but at least from the theoretical point of 
view one would like to understand the choices involved. Furthermore, in some cases 
discussed in the literature, such as models involving potentials with sharp "features" 



r^-[^, the polynomial approximation is by definition unlikely to be sufficient. 



In view of the above it becomes interesting to consider alternative schemes of 
dealing with the flow equations. The purpose of this note is to offer a method of solv- 
ing the hierarchy, which does not set an infinite number of flow parameters to zero. 
Indeed, all the flow parameters are non-zero in this approach, and in consequence 
this method does not restrict the space of generated potentials to polynomials in the 
infiaton field. The hierarchy effectively terminates because flow parameters of higher 
order are expressed algebraically in terms of a finite number lower order ones as a 
consequence of a condition which requires that a flow parameter of some order be 
constant. 
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The immediate question is then how the new set of potentials is related to the set 
of polynomial potentials scanned in Kinney's procedure. It is obvious from the new 
termination condition (introduced in section ^ that the space of potentials accessible 
with the new method is strictly a superset of the old one. One way to see the impact 
of the new scheme would be to repeat Kinney's numerical study and see how the set 
of generated models broadens the regions appearing in [|l| (and in following studies of 
this type). This is clearly interesting to explore. This note however focuses on some 
simple insights which can be gained by analytically solving the new termination 
conditions at low levels, very much in the spirit of Liddle's analytic solution of 



the original scheme [|1^]. Unlike that case however, a complete analytic solution 
is possible only in the two lowest orders: at higher orders one needs to resort to 
numerical methods. The analytic solutions described in section ^ show that the new 
scheme, in accordance with expectations, brings in non-polynomial potentials. In 
particular, the lowest order solution describes power law inflation, which is outside 
the standard scheme since it requires an exponential potential. The next order 
leads to some interesting cases which have appeared in the literature in various 



contexts |23]-[25]. They include so-called "ultra-slow-roll inflation" [26|, which has 
the second Hubble parameter rjH = "i- At higher orders one can write down some 
special solutions analytically, but for practical applications numerical integration is 
required. One can proceed to integrate the flow equations directly, or alternatively 
integrate the termination condition, which is a single ordinary differential equation 
of order M + 1, where M is the order of the flow parameter required to be constant. 
The second approach directly gives H{(j)), the Hubble parameter as a function of 
the infiaton. Both options can be implemented in a straightforward manner and 
from a technical point of view neither requires anything beyond what is used in 
investigations using Kinney's truncation |1|]. 

The inflationary flow equations, as well as their standard truncation are described 
in section |[ The new scheme is presented in section ^ and some analytic solutions 
of the termination conditions are described in section ^. Some closing comments are 
offered in section ||. 



3 



2. Truncating the Flow Equations 



The inflationary flow equations introduced in [|T| assume that inflation is driven by 
a single scalar field described by an effective action of the form 



s = -j d'x^{-{d<j>y + vm . (2.1) 

For spatially homogeneous field configurations Einstein equations reduce to 



-3Hip + p) 



where 



P ■ 

3Mj>H'^ = p 



P = \4>^ + V{4>) 



(2.2) 
(2.3) 



(2.4) 
(2.5) 



Here Mp is the reduced Planck mass (Mp = I/SttG), the dot indicates a time 
derivative and H = a/a. 

It is convenient to write these equations in first order form, treating (p as the 



evolution parameter in place of t. From ( p.2|) - ( p.5|) it follows that 



-2MpH'{(f)) 



(2.6) 



where the prime denotes a derivative with respect to (f). Using this and (|2.5| ) in (|2.3| ) 
gives 

2Mj,H\cj)f = 3Mj,H\(f)) - V{(j)) . (2.7) 



This is the Hamilton- Jacobi form of the field equations p7[]-[p0[]. 

The fundamental indicator of infiation is the first Hubble fiow parameter 

en = 2MU^r , 



:2.8) 



where the prime indicates a derivative with respect to the infiaton field. The basic 
property of en is that 

- = H'il - en) , (2.9) 
a 
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which shows that the Universe is inflating if and only ii en < ^■ 

The number of e- folds at some time t before the end of inflation at time t/ is 
given by 

N = [ ^ Hdt , (2.10) 



so one has dN = —Hdt. This convention defines as the number of e-folds before 
the end of inflation at iV = 0. Thus as time flows forward, decreases. From ( [^.1U| ) 
it follows that 

— = -i- (2 11) 

dN HdcP' ^ ^ 

which can be rewritten using (|2.6|) as 



— = 2M| — — . 2.12 

dN ^ H deb ^ ' 



By direct computation one then finds 

den 



dN 

where 



-2eH{eH - Vh) . (2.13) 



7]H = ^Mj,^ (2.14) 
II 

is the second flow parameter. 

The derivative of rjn involves the third derivative of H, which motivates the 
introduction of another dimensionless flow parameter. Proceeding in this way all 
higher derivatives of the Hubble parameter appear and an infinite hierarchy of dif- 
ferential equations is generated. It can be described compactly by introducing the 
infinite sequence of Hubble flow parameters [|I| defined as 



H 



A. = (2AfJ)'<^^ . > 1 . (2.15) 

SO that Ao = eH and Ai = rjH- The flow equations can now be written as 
(iAo 



dN 
dX 



k 



dN 



2Ao(Ao-Ai) (2.16) 
( - A;Ao + (A; - l)Ai) A^ + A^+i , k>l . (2.17) 
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This is an infinite hierarchy of differential equations for the Hubble flow parameters 
Afc. Solutions of these equations for which Aq < 1 for a sufficiently long time describe 
inflating spacetimes of interest in cosmology. 



It was emphasized by Liddle |jT6[ that the flow equations do not reflect any specific 
choice of potential, since their derivation does not make use of the Hamilton- Jacobi 
equation. This is consistent with their purpose, which is to describe inflationary 
solutions for canonical scalar field theories without prejudice. While the equations 
themselves do not involve a choice of scaler potential, any specific solution of ( p. 17] ) 
corresponds to a specific scalar field theory. This is because once en is found as 
a function of by solving the flow equations one can calculate H{(f)) (up to an 
overall scale) ^. This in turn determines the scalar potential via the Hamilton- Jacobi 
equation (^.7]): 

V{(j)) = mlH{(t)f - 2M^pH\(t)f . (2.18) 

Thus every solution of the flow equations determines the corresponding scalar po- 
tential (up to an overall energy scale). 

The set of all solutions of the flow equations is identical to the set of all solutions 
of the Hamilton- Jacobi equations for all choices of V^(0). In that sense the full space 
of solutions does not reflect any choice of dynamics - it conveniently parameterizes 
the outcome of all the possible choices. Choosing a class of solutions (a subspace of 
all solutions) is however tantamount to a statement of dynamics, and this is what 
practical applications of the flow equations do. 

The procedure introduced by Kinney and elaborated on by many authors, 
involves truncating the infinite hierarchy by setting, for some integer M, 

Afe = , k>M. (2.19) 

This yields a closed set of differential equations for Ao...Aa/. It is important to 
note that truncating the flow equations is not an approximation: solutions to the 
truncated set of equations are exact, but they span a subset of all the solutions to 
the flow equations. Thus truncation is equivalent to restricting the set of all possible 



This also requires relating N and 0, which can be done using dN = 2Ml§Td4>, which follows 
from ( pl0| ) and ( ^ . 
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potentials to some subset. This fact was made explicit by Liddle [|T^ who observed 
that the truncation condition (|2.19|) could (using (|2.15|) ) be written as 



^ - , (2.20) 

which makes it plain that solutions of the flow equations are polynomials of order 
M: 

M 

i/(0) = J]afc0^ (2.21) 

A:=0 



Using this in (|2.18| ) implies that the corresponding scalar potentials are polynomials 



in of order 2M. One can survey a large space of potentials by truncating the flow 
equations at a high level, i.e. by taking M large in ( p.l9| ). 



3. New Solutions of the Flow Equations 



While the class of polynomial potentials appearing in the standard treatments of 
the flow equations may be sufficient for most practical purposes, from a theoretical 
perspective it seems somewhat limited. Indeed, from the point of view of embedding 
inflationary scalar field theories in string theory it seems that this restriction is quite 
severe, since non-polynomial contributions to scalar potentials are quite common in 
that setting |T^, [TB[. It turns out however that a very simple modification of the 
standard truncation of the fiow hierarchy significantly broadens the set of potentials 
covered without introducing any significant complications relative to the standard 
procedure outlined in the previous section. 

The idea is to replace the truncation condition ( |2.19|) by 

Am = A (3.1) 

for some M, where A is a constant. This termination condition closes the fiow 
equations hierarchy at level M by introducing the constant A. The hierarchy closes, 
because only the first M differential equations are non-trivial if (|3.1|) is imposed. The 
equations at levels M and above become algebraic. The termination condition ( |3.1j ) 
does not set the higher order fiow parameters to zero: they are instead expressed in 
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terms of the lower order parameters. For example (|2.17 ) and imply 



Am+i = a(mAo-(M-1)Ai) . (3.2) 

Similar relations can be written down for higher flow parameters which generically 
remain non- vanishing. 

The flow equations can be integrated as before for any choice of M and some 
reasonable set of values of A. The original subset of inflationary model space is the 
case of A = 0, so clearly all the solutions appearing in the old approach are recovered. 

There are in fact two ways to proceed. One option is to integrate the set of 
M + 1 nontrivial flow equations. The alternative is to directly solve the termination 
condition itself. In the case of standard truncation the solution of the truncation 
condition (|2.19|) is ( p.21j ). This way H is obtained directly, without going through 



the flow parameters, in fact circumventing the flow equations themselves. In the 
case of the modified termination condition one can proceed in the same spirit by 
expressing (|3.1|) using ( p.l5|) as 



( TTl\M-l ^M+l TT 



This is a single differential equation of order M + 1 which replaces Liddle's ( |2.2CI| ) 



For M > 1 this equation is nonlinear, and one cannot solve it analytically. It is 
however straightforward to solve numerically. For that purpose it is convenient to 
write it as a system of first order equations as follows. Introducing 

H,^^ fc = 0,...,M (3.4) 

equation ( p.3|) can be rewritten as a system of first order differential equations: 

^ = Hu+, k = 0,...,M-l 

M _ A Hq 



dHjir A 



(3.5) 



Supplementing (|3.5|) with suitable initial conditions one can generalize numerical 
computations of the type pioneered by Kinney [Q to the wider set of solutions de- 
scribed here. To this end one can express initial values for the Hk in terms of often 
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used initial values for the standard flow parameters. Indicating initial values by an 
over-bar one has, from the definition of eu-. 



en 



Hi = ±H,J^. (3.6) 



2Mj> 



The flow equations (or ( p.5| )) determine H up to an overall scale, which can be taken 
as Hq- The choice of sign above reflects the possibility of the inflaton rolling to the 
left or to the right. Similarly for k > 1 one can write 



Hk+i - (^/HoX,^ ^ , (3.7) 



where are initial values of the flow parameters, which can be related directly to 
those used by Kinney |l|]: one has Ai = ^{(Th + ^(^h) where au is the initial value 
of Kinney's an, and Afc are the initial values of Kinney's (for k > 1). One can 
now numerically integrate the equations ( |X3| ) choosing initial values of en, cth, ^^h 
from the same ranges as those used in |^ (and most of the literature devoted to this 
subject) to facilitate comparison. 

The numerical computations following from this prescription will not be pre- 
sented here; instead the following section will describe some analytic considerations 
which give a glimpse of space of potentials defined by the procedure introduced above. 



4. Some Analytic Results 

To understand the difference in the space of potentials scanned by the truncation of 
the flow equations described in the last section it is instructive to look at the simplest 
cases, that is when ( p.l|) is imposed with M = and M = 1, which are very simple 
to solve analytically. 



In the case M = the termination condition ( p.l| ) is the statement that Aq 
is constant. Since Aq is just en-, this is power law inflation (when A < 1). The 



termination condition (|3.1|) becomes 



2Ml{^f = A . (4.1) 

For this to make sense one can only allow non-negative values of A. As discussed in 
the last section, (O) can be regarded (in the spirit of [ JT|) as a first order differential 
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equation for H{(f)). The general solution is 



A 



H{<p) = Aexp{±^l—^<P) . (4.2) 



This involves one integration constant, A. The corresponding potential, obtained 
form ( p.l8| ) is 



V{(P)=A'MUS-X)exp{±^l^^) . (4.3) 

This is the well known example of Lucchin and Matarrese |^ . 

This simplest case already shows the difference between the procedure proposed 
in the previous section and the one normally used in the literature. While the 
standard procedure is equivalent to scanning over the set of polynomial potentials 
of some order, here one obtains a non-polynomial one. Clearly, all the Hubble flow 
parameters are non-zero: they are all given by powers of the constant A. The case of 
standard level truncation is obtained in the limit A ^ 0, which describes de Sitter 
expansion. 

Imposing the termination condition (|3.1| ) with M = 1 is also solvable, and rather 
interesting. The termination condition reads 

H" 

2Ml— = A . (4.4) 
H 

Here there is no restriction on the sign of A, and the character of the solutions of this 
differential equation depend on this sign. Since the differential equation is of second 
order there will be two integration constants. The general solution is^ 



v4cosh(^^0) + 5sinh(y'^0) for A > 

A + B(f) for A = (4.5) 

^cos( J^0) +Bsm{J]^<p) for A < . 



The case A = is of course the result of the standard truncation at this level. 

For a sensible inflationary solution one has to make a choice of integration con- 
stants and restrict the range of (p appropriately so as to ensure that H' does not 



■'Potentials of this type were previously considered in references [ p3| , They can all be char- 
acterized as models with constant rjH- 
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change sign. Rather than discuss this further, this presentation will focus on one 
special case, that of A > with B = 0: 



H{<P) = AcosH^I^<ly) . (4.6) 



2M^ 



The Hubble slow-roll parameters are 




2M|, 



(4.7) 

Vh = X. (4.8) 
The corresponding potential, which follows from ( |2.18| ) reads 



y(0) = ^'M|(A + (3-A)cosh2(^^0)) . (4.9) 

Note that for A = 3 the potential becomes constant: this solution was discussed 
by Tsamis and Woodard ^6) under the name ultra-slow-roll inflation. It was later 
analysed by by Kinney as an example where the spectrum of curvature perturbations 
is exactly scale invariant but where the horizon-crossing formalism fails |2^. The 
solution ( |4.6| ) is valid for a range of A, so taking A close to 3, but not exactly 3, 
should provide interesting examples with almost scale invariant spectra^. Since the 
potential ( [4.9|) in this case is not polynomial, such examples are outside the realm of 
standard truncated flow equation simulations. 

Terminating the hierarchy using (|3.1|) with M > 1 involves solving a nonlinear 
equation of order M + 1. The general solution, depending on M + 1 constants of 
integration appears not to be available analytically, but one special solution is easy 
to write down: it is 

H{<P) = Aexp(±-i=Aw<^) . (4.10) 

ZlVlp 

Both choices of sign are admissible, but one cannot take linear combinations, as the 
equation is not linear. Therefore for higher-level solutions one needs to resort to 
numerical integration, as discussed in the previous section. 



■^A very interesting and nontrivial class of potentials with exactly scale invariant spectra was 
obtained recently by Starobinsky . It appears that these models are not of the type considered 
here. 
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5. Conclusions 



The inflationary flow equations are the basis of a very widely used approach to explor- 
ing the realm of inflationary scalar fleld theories. The standard method of truncating 
the inflnite hierarchy of flow equations restricts the class of scalar potentials to poly- 
nomials in the inflaton fleld. This paper presented a different, more general, way of 
solving the hierarchy. The resulting class of potentials includes those obtained by 
the standard truncation method, but is much broader in that it also includes a wide 
range of non-polynomial potentials. The procedure boils down to solving the system 
of differential equations ( |3.5|) , which is the main result presented here. 

Some insight into the space of potentials accessible using this method can be 
gleaned from analytic solutions to lowest level termination conditions. As discussed 
in section ^, examples found this way have already appeared in the literature in 
various contexts. Here they serve the purpose of illustrating how the new solutions 
of the flow equations extend the range of potentials scanned. 

The approach described here could be particularly useful for exploring models of 
inflation involving potentials with "features" , where the potential is non-polynomial 
in an essential way |[T^-[^. It would also be interesting to understand the impact of 
the present work on numerical flow equation simulations used to analyze data from 
WMAP and other sources g, |0H3 . 



As this paper was being written up the generalization of the flow equations to the 
case of DBI models was introduced by Peiris et al. [|l3l . The new procedure presented 
here can be used without change also in that case. Due to the occurrence of the 
"Lorentz" factor 7 in the Dirac-Born-Infeld action, the hierarchy of flow equations 
requires two truncation conditions, which Peiris et al. solved, demonstrating that 
their method scans the space of V{(f)) and 7(0) which are polynomials in the inflaton 
fleld. Replacing the two truncation conditions used in |[T^ by conditions of the sort 
advocated here works in the same way as the canonical case discussed in this note. 
It would clearly be of interest to investigate how modifying that study along these 
lines affects the results reported there, given the tight observational constraints on 
that class of models. 
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